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ABSTRACT 

We consider new Suzaku data for NGC 3516 taken during 2009 along with other recent X-ray 
observations of the source. The cumulative characteristics of NGC 3516 cannot be explained 
without invoking changes in the line-of-sight absorption. Contrary to many other well-studied 
Seyfert galaxies, NGC 3516 does not show a positive lag of hard X-ray photons relative to soft 
photons over the timescales sampled. In the context of reverberation models for the X-ray 
lags, the lack of such a signal in NGC 3516 is consistent with flux variations being dominated by 
absorption changes. The lack of any reverberation signal in such a highly variable source disfavors 
intrinsic continuum variability in this case. Instead, the colorless flux variations observed at high 
flux states for NGC 3516 are suggested to be a consequence of Compton-thick clumps of gas 
crossing the line-of-sight. 

Subject headings: galaxies: active - galaxies: individual: NGC 3516 - galaxies: Seyfert - X-rays: galaxies 



1. Introduction 

Signatures of ionized gas are commonly ob- 
served in Seyfert galaxies. Well-studied sources 
reveal multiple zones of X-ray absorbing gas cov- 
ering a range of ionization state, column density, 
covering fraction and kinematics (e.g. Netzer et al. 
2003; Behar et al. 2003; Kaspi et al. 2002; Steen- 
brugge et al. 2005; Blustin et al. 2005; McKer- 
nan et al. 2007; Blustin et al. 2007). The de- 
tection of deep K-shell absorption lines from very 



highly ionized species of Fe (e.g. Reeves et al. 2004; 
Kaspi et al. 2002; Turner et al. 2008) has shown 
that X-ray absorbing gas can be spectroscopically 
traced up to equivalent hydrogen column densi- 
ties iVn ^ lO^^'cm^^. The possible existence of 
circumnuclear gas at high column densities is of 
great interest with regard to understanding the 
mass flow around the black hole and energetic 
considerations for the system. UV spectroscopy 
unequivocally shows AGN to possess complex ab- 
sorption systems and so the presence of a complex 
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of X-ray absorbers should be unsurprising. 

The systematic spectral hardening exhibited 
by Seyfert-type AGN as they drop in X-ray flux 
(e.g. Papadakis et al. 2002; Pounds et al. 2004a,b; 
Vaughan & Fabian 2004; Miller et al. 2007, 2008; 
Turner et al. 2008) can be modeled by changes in 
covering fraction of X-ray absorbing gas. There 
are a few cases where changes in individual lines 
trace rapid changes in the X-ray absorber. For 
example, an absorption line from Fe xxv detected 
in NGC 3783 varies on timescales of days (Reeves 
et al. 2004), likely originating ~ 0.1 pc from the 
nucleus. In NGC 1365 variable absorption lines 
have been detected from Fe XXV and Fe XXVI sup- 
porting a picture in which the nucleus suffers vari- 
able obscuration by an absorber whose covering 
fraction changes on short timescales (Risaliti et al. 
2005, 2007). 

In the broad-line Seyfert 1 galaxy (BLSyl) 
NGC 3516 (z=0.008836; Keel 1996), X-ray data 
covering 0.5-10 keV have revealed a strong signa- 
ture from a variable X-ray absorbing outflow. X- 
ray grating data from Chandra HETG and XMM 
RGS show discrete absorption features tracing a 
range of column densities, ionization states and 
velocity components for the gas (e.g. Turner et al. 
2008; Mchdipour ct al. 2010). Here we present 
new Suzaku data from NGC 3516, taken during 
2009 October. We also reconsider recent X-ray 
observations of this source, seeking to reconcile 
the spectral and timing behaviour and thus elu- 
cidate the true nature of the X-ray signatures of 
this AGN. 

2. The Suzaku Observations 

Four co-aligned Suzaku (Mitsuda ct al. 2007) 
telescopes focus X-rays on to CCD cameras com- 
prising the X-ray Imaging Spectrometer (XIS 
Koyama et al. 2007). XIS units 0,2,3 arc front- 
illuminated (FI) while XIS 1 is a back-illuminated 
CCD. XIS 1 has an enhanced soft-band response 
but lower area at 6 kcV than the FI CCDs as 
well as a larger background level at high energies. 
XIS 2 failed on 2006 November 9 and hence was 
not used. Suzaku also carries a non-imaging, col- 
limated Hard X-ray Detector (HXD, Takahashi 
et al. 2007), whose PIN detector provides useful 
AGN data typically over the range 15-70 keV. 

A Suzaku observation of NGC 3516 was made 



2009 Oct (OBSID 704062010) and the data were 
reduced using v6.9 of HEAsoft. We screened the 
XIS and PIN events to exclude data during pas- 
sage through the South Atlantic Anomaly, start- 
ing and ending within 500 seconds of entry or exit. 
Additionally we excluded data with an Earth el- 
evation angle less than 10° and cut-off rigidity 
> 6 GeV. The source was observed at the nom- 
inal center position for the HXD. FI CCDs were 
in 3 X 3 and 5x5 edit-modes, with normal clock- 
ing mode. We selected good events with grades 
0,2,3,4 & 6 and removed hot and flickering pix- 
els using the SISCLEAN script. The spaced-row 
charge injection (SCI) was utilized. The exposure 
times were 222 ks per XIS unit and 178 ks for the 
PIN. XIS products were extracted from circular re- 
gions of 3.5' radius while background spectra were 
extracted from a region of radius 2.5' (offset from 
the source and avoiding the chip corners, where 
there are data from the calibration sources). 

For the 2009 HXD PIN analysis we used the 
model "D" background (released 2008 June 17 ^). 
The time filter resulting from screening the ob- 
servational data was applied to the background 
events model. The ftool HXDPINXBPI was used 
to create a PIN background spectrum from the 
screened background data, hxdpinxbpi takes ac- 
count of the form and flux level of the cosmic 
X-ray background (Boldt & Leiter 1987; Gruber 
et al. 1999) in the Suzaku PIN field of view. We 
used the response file 

ae_hxd_pinxinome3_20090826 . rsp for spectral 
fitting. 

During the 2009 Suzaku observation NGC 3516 
was found to have source count rates 0.507 ±0.001 
(summed XIS 0,3 over 0.75-10 kcV) and 0.049 ± 
0.002 (PIN over 15-50 keV) ct/s. The background 
was 5% of the total XIS count rate in the full band. 
The source comprised 9% of the total counts in the 
PIN band. The source fluxes for recent observa- 
tions of NGC 3516 are shown in Table 1. 

Spectral fits utilized data from XIS and 3, in 
the energy range 0.75 — lOkeV and from the PIN 
over 20 — 50keV. Data in the range 1.75-1.9kcV 
were also excluded from the XIS spectral analysis 
owing to uncertainties in calibration around the 
instrumental Si K edge. We note that data can be 



http;//www. astro. isas.jaxa.jp/suzaku/doc/suzakumemo/ 
suzaJcumemo-2007-01.pdf 
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used down to lower energies, 0.3 keV, for timing 
analysis. XIS 1 was not used owing to the higher 
background level at high energies. XIS data were 
binned at the HWHM resolution for each instru- 
ment, optimal for detection of spectral features, 
while PIN data were binned to be a minimum of 
5a above the background level for the spectral fit- 
ting. In the spectral analysis, the PIN flux was 
increased by a factor 1.18 for the new data, ap- 
propriate for the cross calibration of instruments 
for the 2009 observation ^. The new data were 
considered in conjunction with those from 2005, 
whose screening and data reduction are described 
by Markowitz et al. (2008). 

3. Spectral Fitting Results 

Chandra images of NGC 3516 (George et al. 
2002) have shown an extended soft X-ray emis- 
sion component contributing an observed flux 
-fo.4-2kcV ^ lO^^^'ergcm"^ s~^ within ^ 2 kpc 
(assuming Ho = 73kms~^kpc~^ throughout). 
In addition, there is a known off-nuclear X-ray 
source, CXOU 110648.1-1-723412 (George et al. 
2002) that cannot be separated from the nuclear 
emission using Suzaku. However these two provide 
only a small fraction of the total soft-band flux for 
the epochs considered here. 

For an immediate visual comparison of data 
from recent X-ray observations of NGC 3516, we 
plotted the 2009 XIS spectral data with the 2005 
Suzaku observation and other observations made 
using XMM overlaid (Figure 1). This comparison 
reveals that during 2009 Suzaku found NGC 3516 
in a relatively low flux state. The data show that 
the linc-of-sight absorption is obviously changing 
between epochs for NGC 3516. Although the 
hard X-ray flux levels are very similar in the 2005 
Suzaku and 2001 Apr XMM data, the former shows 
a large amount of extra absorption compared to 
the latter. The lowest level of hard X-ray flux is 
observed in the new data from 2009, although at 
this epoch the absorption appears to be also rel- 
atively low. It is clear that the source cannot be 
simply described by a response in ionization of the 
absorbing gas to changes in the continuum flux. 



^ ftp: / /legacy, gsfc. nasa. gov / suzaku /doc/ xrt / suzakumemo- 
2008-06.pdf 
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Fig. 1.— The recent X-ray states of NGC 3516 
from XMM observations in 2001 and 2006 and 
from Suzaku observations during 2005 and 2009. 
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Fig. 2.— Suzaku XIS and PIN data from 2005 
(blue) and 2009 (red) along with the absorption 
model (black solid line). 
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3.1. Fitting Absorption Models 

From detection of discrete absorption lines and 
from broad spectral curvature, Turner et al. (2008) 
found evidence for a multi-zoned X-ray absorber 
in a joint XMM/ Chandra campaign during 2006. 
Markowitz et al. (2008) also found complex ab- 
sorption to be required to fit the spectrum of 
NGC 3516 observed by Suzaku during 2005 (and in 
that spectral model, additionally required a mod- 
estly broadened component of FeKa emission). 

We fit the Suzaku 2005 and 2009 data over 0.75- 
50.0 keV using xstar v2.2 to parameterize the ion- 
ized gas. The warm absorber table was created 
with effective micro-turbulence velocity dispersion 
a = lOOkms"^ for gas of density n = 10^°cm~'^. 
Here we use the ^ form of ionization parameter, 
^ = L/nr"^, where L is the 1-1000 Rydberg lumi- 
nosity and r the absorber-source distance. 

The spectral model was constructed with a 
powerlaw continuum, adding zones of ionized gas 
until the data were parameterized adequately. The 
photon index, column densities and ionization 
states of the gas were allowed to be free but were 
initially linked to be the same values during 2009 
and 2005. 

We assumed the gas zones to be layered over 
each other, with the highest column density zones 
closest to the nucleus: for this reason, the cov- 
ering fractions do not total 100% (Table 2). The 
model included a full screen of neutral gas fixed at 
the Galactic line-of-sight column, 3.45 x 10^°cm~^, 
covering all components and parameterized using 
the XSPEC model tbabs (Kalberla ct al. 2005). 

The fit yielded an underlying photon index F = 
2.22 ± 0.04 covered by zones of intrinsic absorb- 
ing gas (Table 2) having three different covering 
factors. The outermost layer is Zone 1, compris- 
ing low-ionization gas, as established in NGC 3516 
from Chandra grating spectroscopy (Turner et al. 
2006) and from UV spectroscopy (Kraemer et al. 
2002). HETG spectra required two gas compo- 
nents to accoimt for the soft-band absorption fea- 
tures measured during 2006 (Turner et al. 2008). 
For the joint Suzaku fit we also found that the fit 
was significantly improved (Ax^ = 63) by allowing 
the low-ionization zone to have two components 
(Table 2 Zones la and lb, having the same cover- 
ing fraction). The total column density for Zone 
1 is consistent with the sum of the UV absorbing 



components measured at other epochs (Kraemer 
et al. 2002). 

Additional absorbing zones were found (Ta- 
ble 2) with column densities ~ 2 x lO^^cm"^ and 
~ 2 X lO^'^cm"^, both having log ^ ~ 2. The cover- 
ing fraction of the Compton-thick zone is around 
70% in 2005 and 2009 spectra. We note that this 
hard component can, alternatively, be modeled as 
reflection, as shown by Markowitz et al. (2008). 

XSTAR does not calculate the scattered compo- 
nent of X-ray flux from the absorbing gas, how- 
ever it does calculate the re-emitted line spectrum 
and we included such a component in the model, 
placing it behind the Zone 1 absorber complex. 
With such a placement the re-emission spectrum 
accounts for the strength of Fe Ka emission with- 
out over-predicting the soft line emission (that is 
somewhat suppressed by absorption). The ioniza- 
tion parameter of the emitting gas was left free to 
encompass the possibility that emission arises in 
a gas component other than one of the absorbers 
and indeed, the fltted ionization parameter was 
found to be ^ = 1.53 ± 0.10, different to the ab- 
sorbing gas. The fit gave a reduced Xr = 1-24/358 
degrees of freedom (d.o.f.). Interestingly, the emis- 
sion spectrum fitted here is consistent with the 
ionization-state found for soft-band lines in XMM 
RGS grating data (Turner et al. 2003). 

We stress that the absorption model does not 
over-predict the Fe Ka fiux, or the flux of any 
emission line, even if the gas zones are assumed 
to fully cover the contimium source. The strength 
of the emission from the absorbing gas depends 
on the gas geometry, column density and cover- 
ing/filling factors. Even large column densities 
for the absorber can result in a small predicted 
line equivalent width Miller et al. (2009); Yaqoob 
et al. (2010) when there is expected to be a large 
opacity at the line energy. It is not suprising that 
the absorber emission has not been isolated in 
these data, given the gas geometry suggested in 
this case. 

In this absorption model, the covering frac- 
tions of the gas layers varied significantly between 
epochs (Table 2). Importantly, Zone 2 had a larger 
covering fraction during 2005 (99%) than during 
2009 (86%), despite the hard-band continuum flux 
being a factor of 2 higher (Figure 1). 

The most highly ionized zone detected at any 
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epoch is that traced by absorption Unes from Fc 

XXV and Fe XXVI whose depths indicate an origin 
in gas having A^h ~ 3 x lO^^cm"^ and log^ ~ 4.3 
(Turner et al. 2008) These hncs were detected in 
XMM and Chandra data from 2006, and in Suzaku 
data from 2005, however, they were not detected 
during the Suzaku observation of 2009, with hmits 
on the equivalents widths of (narrow) Fe xxv and 

XXVI absorption found to be < 5 eV and < 20 eV 
respectively. 

The absorption modeling and covering frac- 
tion changes are fundamentally consistent with the 
conclusions of Markowitz et al. (2008) and Turner 
et al. (2008): i.e. the source has a complex ab- 
sorber and changes in covering by a high column 
of gas ~ lO^'^cm^^ of 'intermediate' ionization 
(in terms of the X-ray absorbing zones detected 
in NGC 3516) dominate the spectral variability. 
Small differences between the fit of Markowitz 
et al. (2008) and that presented here arise because 
we are fitting under the joint constraints from 2009 
and 2005 epochs, and also because wo utilize a new 
version of XSTAR and use a slightly different model 
construction. 

Constraining the covering fractions of the gas 
to be the same for 2005 and 2009, allowing only 
the column densities and ionization states of the 
layers to be free for each epoch, does not provide 
a satisfactory joint fit to the two epochs (Xr = 
4.8/355 d.o.f.). It is clear that the source cannot 
be parameterized satisfactorily without invoking 
such a mode of variation. 

Small improvements to the fit can, of course, 
be obtained by allowing the absorber column den- 
sities and ionization-states to also vary between 
epochs. However, as Suzaku cannot distinguish 
which of the layers of gas are varying, we did not 
find it useful to pursue this line of consideration 
further. 

3.2. Fe Ka Line Variability 

To test for variability in the Fe Ka line emission 
we re-fit the 4.0-7.5 keV data using a simple ab- 
sorbed power-law model for the local continuum. 
The line was modeled using a simple Gaussian pro- 
file. Initial fits found the line width to be consis- 
tent with the HEG value tr = 40 eV, at all epochs 
(Turner et al. 2008) and so the fits were rerun with 
the line width fixed at that value. 



The joint 2005/2009 Suzaku fit yielded a peak 
energy for the line E= 6.40 ± 0.06 keV. The nor- 
malization of the FeKa line was 5.75 ± 0.24 x 
10~^ photons cm~^ s""'^ (an equivalent width of 
164 eV against the observed continuum) dur- 
ing 2005 and 4.16to 32 x 10"^ photons cm'^s"^ 
(equivalent width 253 cV) during 2009. Thus com- 
parison between the two mean spectra observed 
using Suzaku reveals strong evidence (> 99.9 % 
confidence) for variability in FeKa line flux over 
timescales of years. 

To investigate the line behavior on shorter 
timescales, the 2001 Apr, 2001 Nov and 2006 
XMM data were re- fit with the new model. Fe Ka 
line fluxes are detailed in Table 3. In addition to 
these fits, the individual OBSIDs comprising the 
XMM data for 2006 were fit individually, yielding 
Fe Ka fluxes sampled on timescales of days: how- 
ever, no significant variability was detectable on 
such a timescale (Table 3). 

4. Timing Analysis 

We examined the hard (4.0 - 7.5 keV) and soft- 
band (0.3-1 keV) light curves from the Suzaku ob- 
servations. As was evident from the spectral anal- 
ysis (Figure 2), the hardness ratio is very different 
for the two epochs (Figure 3). 

The best constraints on the timing properties 
of NGC 3516 are currently offered by the com- 
bined data of 2006. At that epoch the source 
was bright and, taken together, the XMM and 
Chandra exposures provide a long baseline of cov- 
erage (Turner et al. 2008). Wc re- constructed 
the combined XMM/ Chandra light curve in se- 
lected energy bands (Figure 4). The effective areas 
and spectral responses of the XMM and Chandra 
data differ: to create a single light curve, the ex- 
pected XMM and Chandra count rates in narrow 
energy bins were calculated from the mean spec- 
tral model. Next, the time-dependent observed 
XMM count rate in those energy bins was rescaled 
to the expected Chandra count rate using that 
mean conversion between instruments. The light 
curve in broader energy bands was then calculated 
by summing the renormalized values, propagat- 
ing XMM statistical errors, taking account of the 
energy-dependent re-normalization. By scaling to 
the Chandra instrument the uncertainties on the 
Chandra part of the light curve remain at their 
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Poisson values, whereas on the XMM part the un- 
certainties are shghtly worse than Poisson, albeit 
still much smaller than the Chandra uncertainties. 
By re-normalizing in fine energy bins, to first-order 
the efi^ect of spectral variations on the accuracy of 
the cross-instrument combination are minimized. 

For a quantitative assessment of the source vari- 
ability behavior we follow Miller et al. (2010b, a), 
in which a maximum-likelihood method is used to 
measure both the power spectral densities (PSD) 
and the cross-spectral densities of time series cre- 
ated in a number of broad bands of photon energy. 
We assume that the physical process (es) gener- 
ating the count rate variations in the time series 
is Gaussian (equivalent to the Fourier transforms 
of the intrinsic time series, before being obser- 
vationally sampled, having uncorrelated phases). 
The PSD and cross-spectral density are defined 
as the Fourier transforms of the autocorrelation 
and cross-correlation functions. To find the best- 
fitting PSD and cross-spectra we compute the ex- 
pected covariance matrix between time samples 
for an initial set of PSD and cross-spectral values 
in bins of frequency and iterate to a maximum like- 
lihood solution using the method of Bond, Jaffe, & 
Knox (1998). The covariance matrix includes the 
effects of shot noise. The time lags between bands 
are evaluated as parameters in the cross-spectral 
density (Miller et al. 2010b,a). 

We applied these methods to the 2006 com- 
bined XMM/ Chandra time series; above lO""' Hz 
the XMM data were used alone, owing to the 
low signal-to-noise ratio for Chandra data sam- 
pled at those frequencies. This new evaluation of 
the PSD (Figure 6) confirms previous conclusions 
(Edelson & Nandra 1999) that the PSD shows a 
steep drop to high frequencies. In contrast to the 
long XMM/ Chandra campaign of 2006, the Suzaku 
data yielded a relatively short baseline, relatively 
low amplitude of variability and did not yield con- 
clusive results for the lag spectra. Intriguingly the 
lag spectra for 2006 data from NGC3516 showed 
no positive time lags (Figure 5). At the lowest fre- 
quencies sampled there is weak evidence for neg- 
ative time lags (soft band lagging the hard band) 
but as the timescales concerned are approaching 
the longest periods sampled by these observations, 
caution is required in their interpretation, and in- 
dependent observations are required to confirm 
their reality. 



, („ V>", '". , " ""i* ' ' i, 

5x10' 10' 1.5x105 2x10= ° 105 2^,05 3,^05 4,,os 

TIME s RATE count/s 

Fig. 3. — Suzaku XISO+3 time series in 2048s bins, 
from the 2005 and 2009 Suzaku observations. Red 
points represent 4-7.5 keV and black the 0.3-1 keV 
band. 
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Fig. 4. — The light curve in 2048s bins from 
2006 XMM pn and Chandra HETG data with soft 
(black: 0.3-1.0 keV) and hard bands (red: 4.0- 
7.5 keV) overlaid. The XMMlight curve rates have 
been renormalized to the Chandra count rate, and 
HETG points can be recognized most easily by the 
relatively large uncertainties on the count rate. 
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5. Discussion 



10-' 



10 = 



Fig. 5. — Lag spectra for 2006 XMM data, calcu- 
lated for 4-7.5 keV versus 0.3-1 keV. On this plot, 
positive values of r would represent a lag of hard 
photons relative to the soft. 




Fig. 6.— PSD from 2006 data over 0.3-1.0 keV. 
The open points denote the combined XMM and 
Chandra data using 2048 s bins, the solid points 
denote the XMM data alone, sampled using 96 s 
bins. 



5.1. Spectroscopy Results 

Comparison of recent X-ray spectra of NGC 3516 
shows that the source has a relatively hard spec- 
trum at some epochs of high X-ray flux. A re- 
view of historical X-ray observations of NGC 3516 
confirms that there have been numerous epochs 
during which the source has shown heavy absorp- 
tion by low-ionization gas, often at high levels 
of continuum flux (Kruper et al. 1990; Ghosh & 
Soundararajaperumal 1991; Kolman et al. 1993; 
Guainazzi et al. 2001). The covering changes for 
the absorbing gas are not well-correlated with 
the observed X-ray continuum flux in NGC 3516, 
implying that, in addition to the flux variations 
that are a direct result of obscuration changes, 
one must invoke either intrinsic continuum varia- 
tions or 'colorless' variability originating from the 
passage of Compton-thick blobs of gas across the 
line-of-sight. Reverberation results favor the lat- 
ter possibility (see Section 5.2) and in that case 
the additional component of the model forms a 
natural extension to the absorber spectral model. 

The light-bending model of Miniutti & Fabian 
(2004), which attempts to account for flux and 
spectral variability in AGN by varying the dis- 
tance of the continuum source perpendicularly 
from the accretion disk, cannot naturally account 
for behavior where the source spectrum is found 
to be hard at a high flux state. Adherence to 
this model would require an additional variabil- 
ity mechanism to account for all of the X-ray 
states observed in NGC 3516. Invocation of an- 
other source of continuum variability would re- 
quire an explanation of the lack of positive lags 
between hard and soft bands, which are observed 
in other AGN. 

Other models that have been suggested to ex- 
plain the spectral variability of NGC 3516 are also 
not favored by this comparison of spectral states. 
Netzer et al. (2002) found ASCA observations 
of NGC 3516 to be consistent with changes in 
the ionization state of the X-ray absorber in re- 
sponse to variations in the continuum flux level. 
Mehdipour et al. (2010) also favor the 'responsive 
gas' scenario (although based upon a column den- 
sity an order of magnitude lower than that fitted 
by Turner et al. 2008, owing to their different mod- 
eling of the X-ray continuum and, in particular. 



7 



the inclusion of a black body component that we 
do not require). However, a responsive-gas model 
cannot explain the behavior of NGC 3516 across 
all the X-ray observations; variable-covering ab- 
sorption is required to explain the overall behavior 
of the source. In other work, Bhayani & Nandra 
(2010) have claimed that difFcrcnccs in the RMS 
variability between the 4.6-5.0 keV band and 5.0- 
5.4 and 5.4-5.8 keV bands for the 2001 XMM data 
disfavor an absorption model. However, the pre- 
dicted measure of RMS variability in a given band 
depends on the column and ionization-state of the 
absorbing layer that dominates the variations in 
that band. As the absorption model matches the 
X-ray spectrum of NGC 3516 at all flux levels and 
time slices, it follows de facto that the model re- 
produces also the RMS spectral variability. 

Of course the partial-covering model does not 
rule out a response of the absorbing gas to changes 
in the incident flux: it is inevitable (and therefore 
predicted) that covering changes for inner zones 
of a complex absorber modify the flux and spec- 
tral shape of the continuum reaching zones further 
out, possibly leading to measurable changes in the 
ionization balance of the gas. 

As discussed for the 2006 XMM/ Chandra data 
(Turner et al. 2008), numerous absorption lines 
are detected in the HETG and RGS gratings, pro- 
viding a firm basis for the absorption model in 
NGC 3516. The combined signature of multi- 
zoned X-ray absorbers produces curvature in the 
observed spectral shape, while spectral variability 
can be explained by changes in covering fraction 
of the layers. Although current X-ray data do not 
allow the gas covering fraction to be measured di- 
rectly from the absorption lines, this quantity can 
be constrained from the broad X-ray spectral cur- 
vature. 

This X-ray absorption model comprises a com- 
pelling extension of the well-studied UV ab- 
sorption systems. NGC 3516 shows some of the 
strongest and most variable intrinsic blueshifted 
UV absorption lines of any Seyfert 1 galaxy (Kol- 
man et al. 1993; Goad et al. 1999). The C iv 
absorption is particularly deep and broad and has 
led some to compare the systems in NGC 3516 
to those in broad absorption line QSOs (Ulrich 
1988). UV absorption features arise in multi- 
ple zones having covering fractions < 1 (Kraemer 
et al. 2002). Kraemer et al. (2002) resolve sev- 



eral kinematic components of absorption in the 
UV lines. Those authors also isolated zones of 
gas having signatures spanning both the UV and 
X-ray regimes, confirming an association between 
the two bands that had been previously suggested 
(Guainazzi et al. 2001). Estimates placed the C 
IV origin at ^ lO-'^^cm (Koratkar ct al. 1996) with 
other lines arising out to ~ lO^^cm from the cen- 
tral source (Kraemer et al. 2002). 

Similarly to Seyfert type 2 AGN, NGC 3516 
shows a bipolar morphology for the naxrow-line- 
region gas (Poggc 1989). In NGC 4151, STIS 
data show evidence; for fast-moving clouds existing 
over a wide angle and originating close to the nu- 
cleus (Hutchings et al. 1999; Crenshaw & Kraemer 
2007), possibly collimating the nuclear radiation 
that excites the NLR gas (Kraemer et al. 2008). 
A similar model may apply to NGC 3516, where 
a toroidal distribution of outflowing clouds might 
coUimate the nuclear radiation. Such a toroidal 
gas distribution could provide X-ray absorption 
and emission features from different zones within 
the flow. If our view of NGC 3516 intercepts the 
edge of the collimating flow we might naturally ex- 
pect the source behavior to be more variable than 
sources viewed at face-on or edge-on orientations. 

In the context of the partial-covering model, 
variations of the Fe Ka flux may arise if the line- 
emitting region suffers variable obscuration. How- 
ever, the line flux also depends on the illuminating 
flux reaching the emitting gas and so the line be- 
havior might be expected to appear complex in a 
multi-layer absorber model. 

HEG data find the Fe Ka emitter to have a bulk 

velocity consistent with zero (Turner et al. 2008). 
The HEG constraint on the line width is 25 < cr < 
50 eV, equivalent to FWHM velocity broadening in 
the range 3000—6000 km/s. For gas in a Keplerian 
orbit the line width corresponds to an origin (4- 
17) sin^i light-days from the nucleus (where i is the 
angle of inclination of the orbital plane, assuming 
a black hole mass 2.95 x 1O''M0, Nikolajuk et al. 
2006). 

5.2. Timing Results 

There has been mounting support for the idea 
that large amplitude and rapid X-ray flux vari- 
ations arc caused by changes in the linc-of-sight 
absorption (Boiler et al. 1997; Guainazzi et al. 
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1998; Tanaka et al. 2004). In such a picture, the 
low/hard states of AGN may be attributable to 
increases in X-ray obscuration, as suggested for 
NGC3516 based on BeppoSAX and ASCA data 
(Nogami et al. 2004). Deep dips in X-ray flux 
have been traced by the light curves of MCG-6- 
30-15 (McKernan & Yaqoob 1998) and NGC3516 
(Turner et al. 2008), where data allow us to track 
the apparent ingress through to the egress of an 
occultation event. Interestingly, the detailed dip 
profiles match the profile of dip events in some 
X-ray binaries, e.g. GRO 1655-40 (Tanaka et al. 
2003), where in that case the dip was attributed to 
an increase in covering fraction of a complex X-ray 
absorber. Spectral variability in Cir X-1 has also 
been attributed to changes in covering of a com- 
plex absorber (Brandt et al. 1996). The deep dips 
in NGC 3516 and MCG-6-30-15 are similarly con- 
sistent with occultation events. The large num- 
ber of similarities in observed timing and spectral 
properties of AGN and stellar binary systems sug- 
gests that complex and variable absorption domi- 
nates the observed X-ray properties for accreting 
systems over stellar to galaxy size scales. This, in 
turn, suggests that the often-outflowing absorber 
is closely related to the accretion process, and 
therefore of the highest importance in understand- 
ing the fundamental mechanisms in these systems. 

Close study of key sources is required to de- 
couple intrinsic continuum variations from ab- 
sorption effects. The 2006 XMM/ Chandra cam- 
paign (Turner et al. 2008) showed that in the 
0.5-10 keV band, there appear to be two char- 
acteristically different modes of X-ray variabil- 
ity: one which results in opacity changes and 
thus changes in hardness ratio, the other, partic- 
ularly at high fluxes, which results in no change 
in hardness ratio. In the context of the com- 
plex absorption model, the latter mode may ap- 
ply when the continuum is uncovered by the fitted 
absorbers: in the 2006 data this occurred above 
a flux F2-10 ~ 3 X 10~^^ergcm~2s~-^. Variability 
above this threshold was attributed either to in- 
trinsic variations in the continuum or to changes in 
covering by Compton-thick 'bricks' (Turner et al. 
2008). 

Recent long X-ray exposures on two bright 
AGN have allowed signiflcant progress with re- 
gard to our general understanding. Suzaku data 
from NGC 4051 show that the frequency- and 



energy-dependence of the observed time lags can 
be explained by the effects of reverberation from 
a thick shell of material extending to ~ 600 grav- 
itational radii from the continuum source (Miller 
et al. 2010b). That work was the first interpre- 
tation of the previously-known positive AGN time 
lags, that dominate the lag spectra, as arising from 
reverberation within the X-ray band. 

A second important result arises in XMM data 
from IH 0707-495. The source shows lags that os- 
cillate between positive and negative values as a 
function of temporal frequency of variation. How- 
ever, in IH 0707-495 it has also becin shown (Miller 
et al. 2010a) that a simple reverberation model 
can explain the full set of (positive and negative) 
lags, requiring no additional physical processes to 
explain the observed lag sign changes. 

Many AGN show a lag of hard X-ray photons 
relative to soft, with the lag increasing to decreas- 
ing temporal frequency (e.g. McHardy et al. 2004; 
Papadakis et al. 2001; Markowitz 2005; Arevalo 
et al. 2006; Dasgupta & Rao 2006; Markowitz et al. 
2007; Sriram et al. 2009). However, NGC 3516 
shows no such positive lag (Figure 5). In the 
context of our reverberation model for AGN 
lags (Miller et al. 2010b, a), the lack of a lag in 
NGC 3516 is interpreted as an absence of a rever- 
beration signal. 

These imusual lag results may be associated 
with the special sight-line at which we observe 
NGC 3516. The results may be explained if the 
variations observed during 2006 are attributable 
predominantly to covering changes of the absorber 
along our line of sight (i.e. the continuum being 
intrinsically quasi-constant dining those observa- 
tions), as we would not then expect to see sys- 
tematic reverberation, which require all lines of 
sight to vary coherently. Even if the reprocessor 
subtended a large solid angle to the continuum, 
there would be no reverberation for the case where 
there is no intrinsic variation in the continuum. In 
this case, the PIN-band fluctuations might be at- 
tributed to variable obscuration by material hav- 
ing TVh > lO^^cm-^. 

The tentative negative lag observed at low fre- 
quencies in NGC 3516 might be explained by the 
presence of a reprocessed signal in the soft-band, 
rather than in the hard band. Such a soft sig- 
nal would be expected if we observed the emis- 
sion spectrum from Compton-thin, ionized layers 
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of the absorber. This suggestion is supported by 
the presence of significant soft-band hne emission 
at the lowest X-ray flux states. 

6. Conclusions 

Comparison of 2009 Suzaku data with historic 
X-ray observations of NGC3516 sliows that rela- 
tively hard X-ray spectral forms are exhibited at 
some epochs of high X-ray flux. Variable X-ray ab- 
sorption is necessary to explain the range of states 
observed: simple changes in gas ionization state 
are insufficient to explain the range of properties 
exhibited by NGC 3516. 

While much of the variability below 10 keV may 

be explained by changes in the modeled X-ray ab- 
sorber complex, there is a colorless component of 
variability evident during high flux states. In past 
work, the colorless flux variations have been sug- 
gested to originate as either intrinsic continuum 
fluctuations, or to be a consequence of the pas- 
sage of Compton-thick clumps of gas across the 
line-of-sight. 

NGC 3516 provides a rare example to date of 
a source showing no significant lag of hard X- 
ray photons relative to soft photons over the 
timescales studied. This result may be inter- 
preted as an absence of any hard reverberation 
signal. If this interpretation is correct, despite a 
large amount of reprocessing material apparently 
being present around the nucleus, then it seems 
likely that the continuum source is not intrinsi- 
cally strongly variable on the timescales studied. 
We suggest that even the high-state fiux varia- 
tions may be attributed to Compton-thick clumps 
crossing the line-of-sight. The presence of such 
high-opacity clumps is a natural extension of the 
complex absorption model that can describe the 
X-ray spectrum. 

The unusual properties of NGC 3516 may be 
a consequence of us viewing the source at a rare 
orientation, where our line-of-sight intercepts the 
edge of a significant absorber structure that col- 
limates the nuclear radiation leading to the ob- 
served biconical structure for the optical narrow- 
line-region gas. 

TJT acknowledges NASA grant NNX08AL50G. 
This research has made use of data obtained from 
the High Energy Astrophysics Science Archive Re- 



search Center (HEASARC), provided by NASA's 
Goddard Space Flight Center. 

REFERENCES 

Arevalo, P., Papadakis, I. E., Uttley, P., McHardy, 
I. M., & Brinkmann, W. 2006, MNRAS, 372, 
401 

Behar, E., Rasmussen, A. P., Blustin, A. J., Sako, 
M., Kahn, S. M., Kaastra, J. S., Branduardi- 
Raymont, G., & Steenbrugge, K. C. 2003, ApJ, 
598, 232 

Bhayani, S. & Nandra, K. 2010, MNRAS, 408, 
1020 

Blustin, A. J., Kriss, G. A., Holczer, T., Be- 
har, E., Kaastra, J. S., Page, M. J., Kaspi, 
S., Branduardi-Raymont, G., & Steenbrugge, 
K. C. 2007, A&A, 466, 107 

Blustin, A. J., Page, M. J., Puerst, S. V., 
Branduardi-Raymont, G., & Ashton, C. E. 
2005, A&A, 431, 111 

Boldt, E. & Leiter, D. 1987, ApJ, 322, LI 

Boiler, T., Brandt, W. N., Fabian, A. C, & Fink, 
H. H. 1997, MNRAS, 289, 393 

Bond, J. R., Jaffe, A. H., & Knox, L. 1998, 
Phys. Rev. D, 57, 2117 

Brandt, W. N., Fabian, A. C, Dotani, T., Nagase, 
F., Inoue, H., Kotani, T., & Segawa, Y. 1996, 
MNRAS, 283, 1071 

Crenshaw, D. M. & Kraemer, S. B. 2007, ApJ, 

659, 250 

Dasgupta, S. & Rao, A. R. 2006, ApJ, 651, L13 

Edelson, R. & Nandra, K. 1999, ApJ, 514, 682 

George, I. M., Turner, T. J., Netzer, H., Krae- 
mer, S. B., Ruiz, J., Chelouche, D., Crenshaw, 
D. M., Yaqoob, T., Nandra, K., & Mushotzky, 
R. F. 2002, ApJ, 571, 265 

Ghosh, K. K. & Soundararajaperumal, S. 1991, 

ApJ, 383, 574 

Goad, M. R., Koratkar, A. P., Kim-Quijano, J., 
Korista, K. T., O'Brien, P. T., & Axon, D. J. 
1999, ApJ, 524, 707 



10 



Grubcr, D. E., Matteson, J. L., Peterson, L. E., & 
Jung, G. V. 1999, ApJ, 520, 124 

Guainazzi, M., Coniastri, A., Stirpe, G. M., 
Brandt, W. N., Fiore, F., Leighly, K. M., Matt, 

G. , Molendi, S., Puchnarewicz, E. M., Piro, L., 
k Siemiginowska, A. 1998, A&A, 339, 327 

Guainazzi, M., Marshall, W., & Parmar, A. N. 

2001, MNRAS, 323, 75 

Hutchings, J. B., Crenshaw, D. M., Banks, A. C., 
Gull, T. R., Kraemer, S. B., Nelson, C. H., 
Wcistrop, D., Kaiser, M. E., & Joseph, C. L. 

1999, AJ, 118, 2101 

Kalberla, P. M. W., Burton, W. B., Hartmann, 
D., Arnal, E. M., Bajaja, E., Morras, R., & 
Pocppcl, W. G. L. 2005, VizieR Online Data 
Catalog, 8076, 

Kaspi, S., et al. 2002, ApJ, 574, 643 

Keel, W. C. 1996, AJ, 111, 696 

Kolman, M., Halpern, J. P., Martin, C, Awaki, 

H. , & Koyama, K. 1993, ApJ, 403, 592 

Koratkar, A., et al. 1996, ApJ, 470, 378 

Koyama, K., et al. 2007, PASJ, 59, 23 

Kraemer, S. B., Crenshaw, D. M., George, I. M., 
Netzer, H., Turner, T. J., & Gabel, J. R. 2002, 
ApJ, 577, 98 

Kraemer, S. B., Schmitt, H. R., & Crenshaw, 
D. M. 2008, ApJ, 679, 1128 

Kruper, J. S., Canizares, C. R., & Urry, C. M. 
1990, ApJS, 74, 347 

Markowitz, A. 2005, ApJ, 635, 180 

Markowitz, A., Papadakis, I., Arevalo, P., Turner, 
T. J., Miller, L., & Reeves, J. N. 2007, ApJ, 
656, 116 

Markowitz, A., Reeves, J. N., Miniutti, G., Ser- 
lemitsos, P., Kunieda, H., Yaqoob, T., Fabian, 
A. C, Fukazawa, Y., Mushotzky, R., Okajima, 
T., Gallo, L. C, Awaki, H., & Griffiths, R. E. 
2008, PASJ, 60, 277 

McHardy, I. M., Papadakis, I. E., Uttley, P., Page, 
M. J., & Mason, K. O. 2004, MNRAS, 348, 783 



McKernan, B. & Yaqoob, T. 1998, ApJ, 501, L29 

McKernan, B., Yaqoob, T., & Reynolds, C. S. 
2007, MNRAS, 379, 1359 

Mehdipour, M., Branduardi-Raymont, G., & 
Page, M. J. 2010, A&A, 514, AlOO 

Miller, L., Turner, T. J., & Reeves, J. N. 2008, 
A&A, 483, 437 

— . 2009, MNRAS, 399, L69 

Miller, L., Turner, T. J., Reeves, J. N., & Braito, 
V. 2010a, MNRAS, 408, 1928 

Miller, L., Turner, T. J., Reeves, J. N., George, 
I. M., Kraemer, S. B., & Wingert, B. 2007, 

A&A, 463, 131 

Miller, L., Turner, T. J., Reeves, J. N., Lobban, 
A., Kraemer, S. B., & Crenshaw, D. M. 2010b, 
MNRAS, 403, 196 

Miniutti, G. & Fabian, A. C. 2004, MNRAS, 349, 
1435 

Mitsuda, K., et al. 2007, PASJ, 59, 1 

Netzer, H., Chelouche, D., George, I. M., Turner, 
T. J., Crenshaw, D. M., Kraemer, S. B., & Nan- 
dra, K. 2002, ApJ, 571, 256 

Netzer, H. et al. 2003, ApJ, 599, 933 

Nikolajuk, M., Czerny, B., Ziolkowski, J., & 
Gierlihski, M. 2006, MNRAS, 370, 1534 

Nogami, K., Negoro, H., Hong, S., & Mihara, T. 
2004, Nuclear Physics B Proceedings Supple- 
ments, 132, 209 

Papadakis, I. E., Nandra, K., & Kazanas, D. 2001, 
ApJ, 554, L133 

Papadakis, I. E., Petrucci, P. O., Maraschi, L., 
McHardy, I. M., Uttley, P., & Haardt, F. 2002, 
ApJ, 573, 92 

Pogge, R. W. 1989, ApJ, 345, 730 

Pounds, K. A., Reeves, J. N., Page, K. L., & 
O'Brien, P. T. 2004a, ApJ, 605, 670 

— . 2004b, ApJ, 616, 696 



11 



Reeves, J. N., Nandra, K., George, I. M., Pounds, 
K. A., Turner, T. J., & Yaqoob, T. 2004, ApJ, 
602, 648 

Risaliti, G., Bianchi, S., Matt, G., Baldi, A., Elvis, 
M., Fabbiano, G., & Zezas, A. 2005, ApJ, 630, 
L129 

Risaliti, G., Elvis, M., Fabbiano, G., Baldi, A., 
Zezas, A., & Salvati, M. 2007, ApJ, 659, Llll 

Sriram, K., Agrawal, V. K., & Rao, A. R. 2009, 
ApJ, 700, 1042 

Steenbrugge, K. C. et al. 2005, A&A, 434, 569 

Takahashi, T. et al. 2007, PASJ, 59, 35 

Tanaka, Y., Boiler, T., Gallo, L., Keil, R., & Ueda, 
Y. 2004, PASJ, 56, L9 

Tanaka, Y., Ueda, Y., & Boiler, T. 2003, MNRAS, 
338, LI 

Turner, T. J., Kraemer, S. B., Mushotzky, R. F., 
George, I. M., & Gabel, J. R. 2003, ApJ, 594, 

128 

Turner, T. J., Miller, L., George, I. M., & Reeves, 
J. N. 2006, A&A, 445, 59 

Turner, T. J., Reeves, J. N., Kraemer, S. B., & 
Miller, L. 2008, A&A, 483, 161 

Ulrich, M. 1988, MNRAS, 230, 121 

Vaughan, S. & Fabian, A. C. 2004, MNRAS, 348, 
1415 

Yaqoob, T., Murphy, K. D., Miller, L., & Turner, 
T. J. 2010, MNRAS, 401, 411 



This 2-column preprint was prepared with the AAS lAlJ^X 
macros v5.0. 



12 



Table 1 
Observed X-ray Flux 



Observation 


Flux (0.5-2) keV 


Flux (2-10) keV 


Flux (10-50) keV 


2001-04 XMM 


0.49 


2.27 




2001-11 XMM 


0.29 


1.63 




2005-10 Suzaku 


0.13 


2.30 


8.11 


2006-10 XMM/Mean 


1.75 


4.36 




2006-10 XMM/0401 


2.21 


5.08 




2006-10 XMM/0501 


2.02 


4.49 




2006-10 XMM/0601 


1.04 


3.62 




2006-10 XMM/1001 


1.95 


4.44 




2009-10 Suzaku 


0.30 


1.30 


3.90 



Note. — The observed flux in erg cm ^ 



Table 2 
Partial Covering Model 



zone 


Log{ 


JVh ' 


C (2005)^ 


C (2009)^ 


la 
lb 

2 

3 


0.17±0.28 

^■oJ^-o.io 

2.oot?,;°i 
L'.o.T 1 o.or, 


'J-4S_o,16 
„+0.64 
J^-°'-0.17 

20(1 (i 


99 ± 1% 

99 ^ I'/f, 
70 1 2 'A 


97 ± 3% 

86 + 1% 



Note. — A column of neutral gas covered all components, 
representing the Galactic absorption, see text for details. An 

additional column of ionized gas covered all components dur- 
ing 2005 only, with TVh ^ 3 X lO^'^cm"^, log ^ ^ 4.3, see 
text for additional details. Errors arc calculated at 90% con- 
fidence. The remaining flux not accounted for in the Table, 
travels along a line-of-sight that is not covered by the intrinsic 
absorbers. 

^Column density in units of lO^^atom cm"'^ 
^Percentage Covering 



Table 3 
FeK Emission Line 



Observation Normalization^ 



2001-04 XMM 3.33^0 34 

2001-11 XMM 4.02+° 2^ 

2005- 10 S-uzafcM 5.75 ±0.24 

2006- 10 XMM/Mcan 3.94 ± 0.17 
2006-10 XMM/0401 3.51 + o;?3 
2006-10 XMM/0501 i-00Vi%\ 
2006-10 XMM/0601 4.25+°-g3 
2006-10 XMM/1001 3.96l° 
2009-10 Suzaku 4.16toj2 



Note. — For these fits the line energy 
was fixed at 6.4 kcV and the line width at 
a = 40 eV, as determined from HEG data 
(Turner et al. 2008) Errors were calculated 
at 90% confidence 

^Normalization of a Gaussian 
model fit to the Fe Ka line, in units 
10~^ photons cm ~^s~^. 



